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The influence of random interlayer exchange on the phase states of the simplest magnetic het- 
erostructure consisting of two ferromagnetic Ising layers with large interaction radius is studied. 
It is shown that such system can exist in three magnetic phases: ferromagnetic, antiferromagnetic 
and ferrimagnetic. The possible phase diagrams and temperature dependencies of thermodynamic 
parameters are described. The regions of existence of the magnetic phases in external magnetic field 
are determined at zero temperature. 
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Thin films of layered magnets and artificial het- 
erostructures composed of alternating magnetic and non- 
magnetic layers can have a variety of stable magnetic 
states and switching between them can be achieved by 
different regimes of magnetic field variationii2i£. This 
opens vast possibilities for numerous technical applica- 
tions of such structures'^ and makes theoretical de- 
scription of their magnetic states and properties of these 
states very actual. In particular, such description for 
finite number of layers can be achieved in some vari- 
ants of mean-field approximations 6 ' 7 or using numerical 
methods^, yet the study of ideal homogeneous systems 
can be insufficient for the description of experimental 
data. Indeed, impurities and defects in interlayer space 
modify generally the exchange constants, which could 
essentially transform magnetic state of layered struc- 
ture. Fluctuations of nonmagnetic spacer thickness also 
give rise to the similar effect of local interlayer exchange 
modificationiSiii. 

Here we consider the influence of random interlayer 
exchange on the phase states of the simplest magnetic 
heterostructure consisting of two ferromagnetic Ising lay- 
ers. It can be realized as thin two-layer crystal film hav- 
ing controlled concentration of impurities in the inter- 
layer space or as two monoatomic layers deposited on 
the surfaces of nonmagnetic spacer with random thick- 
ness fluctuations. We consider the case when radius of 
intralayer ferromagnetic interaction is much larger than 
lattice parameter. Then the model can be treated in the 
mean-field approximation (except of narrow vicinity of 
transition point), i. e. the intralayer interaction radius 
can be put infinite. In such approximation it is quite easy 
to obtain the expression for averaged thermodynamic po- 
tential depending on the layers' magnetizations, establish 
possible magnetic phases of such system and describe its 
thermodynamic properties in these phases. 



I. THERMODYNAMIC POTENTIAL OF 
RANDOM MAGNETIC BILAYER 

The Hamiltonian of the model has the form 
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Here S n ,i = ±1 are Ising spins, index n numerates 
lattice cites in layers from 1 to N, i = 1,2 is the layer 
number, J > is the constant of intralayer ferromag- 
netic exchange, J n is random interlayer exchange, H is 
external magnetic field. We assume that all J n have the 
common distribution function 



W(J) = (l-p)s(j-J + ^+pS(j-jJ). (2) 



It describes the presence of impurities or defects of one 
sort with concentration p, which change the interlayer 
interaction constant on sites from J + in pure system to 
J_ . We can put J+ > J_ without loss of generality. 

According to Ref. [12j the averaged over J„ inequilib- 
rium potential (per one site in the layer) corresponding 
to Eq. (JU is 



F(m) = Jm 2 /2 
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Here m\ and rri2 are the layers' magnetizations, T - 
temperature, Tr designates the sum over spin configu- 
rations and the angular brackets mean the average over 
random J. Minimizing the potential in Eq. @ over mi 
and m,2 one can find their equilibrium values as functions 
of T, p, H and determine all thermodynamic propeties 
of the model. Summing over spins in Eq. and using 
Eq. we get the explicit form of F(m), 
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F(m) = -Tln2 
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Here to = (mi + TO2) /2 is full magnetization and / = 
(mi — TO2) /2 is antifcrromagnetic order parameter. In 
Eq. Q the unit system is used in which J = 1. 



II. PHASE TRANSITIONS IN ZERO FIELD. 

The potential in Eq. (0J describes generally the com- 
petition between ferromagnetic and antiferromagnetic or- 
dering of the layers. Equating to zero the coefficients at 
m 2 and I 2 in its expansion we find the equations for the 
temperatures of the ferromagnetic transition (T+) and 
the antifcrromagnetic one (T_) correspondingly 



T+ = 1 + (1 -p)tanh(^ 
T_ = 1- (l-p)tanhfi- 



-ptanh ^^r- 
- p tanh ( if- 
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The second order phase transition from paramagnetic 
phase to = 0, / = takes place into the phase which 
have the higher transition temperature. Hence we have 
transition into the ferromagnetic (F) phase for J± > 
and transition into the antiferromagnetic (A) phase for 
J ± < at all p. In both these cases having trivial thermo- 
dynamics with just one equilibrium phase, the disorder 
can induce the appearance of other magnetic phases only 
as metastable ones. This can only cause some peculari- 
ties in the kinetics of the system. 

Naturally, the real competition of two types of equi- 
librium ordering is present when J + > 0, J_ < 0. Then 
there is the critical concentration p c below which ferro- 
magnetic transition takes place while above it antiferro- 
magnetic transition appears. From Eqs. JSJ we have 



tanh( J+) 
P c ~ tanh(./ + )+tanh|J_| ' 

T+ (p c ) = T_ (p c ) = 1. 

Below we focus mainly on this case of competing ran- 
dom interactions. 

Numerical analysis of the F(m) minima shows that on 
the (p, T) - phase diagram the point (p c , 1) is multiphase 
one, in its vicinity three or four magnetic phases can ex- 
ist. The fourth phase is ferrimagnetic (FA) one where 
both to and I are nonzero, so layers have different mag- 
netization modules. The appearance of this phase is the 
main qualitative effect of random interlayer exchange, it 
does not appear without such randomness. 



To elucidate the macroscopic mechanism underlying 
the appearance of such states in the presence of just small 
fraction of negative interlayer bonds one can consider the 
possible ground states describing the properties of the 
system at T = 0. The simple analytical calculations can 
show how spins may arrange to provide the states with 
toi 7^ 771 2 minimizing the energy in Eq. Indeed, from 
Eq. we have for the energy density in some specific 
spin and interlayer bonds configuration 



E = 



/2 + |J_|E_ - J+E H 



Here 
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with the sums running over the sites with negative and 
positive interlayer bonds in £_ and E + correspondingly. 
When negative bonds J_ are absent (p = 0) the lowest 
energy has ferromagnetic spin state with all spins, say, 
up. Introducing negative bonds with some small concen- 
tration could make some spins coupled by such bonds to 
point down in the ground state. Consider the state with 
Ni reversed spins in first layer and -/V2 reversed spins in 
the second one belonging to the different pairs. Then we 
have 

£_ = P — 277l — 2772, E + = I — p, 

rrn = 1 - 2t7,, 77, = Ni/N, i = 1,2, 

with 77; > and rt\ + 772 < p. 

Substituting these expressions in E one can see that 
the resulting E(n\,n2) has no minima inside the region 
where 77; are defined. So the minima are present at the 
boundaries of this region. It is easy to verify that there 
could be two absolute minima, ferromagnetic one (n\ = 
772 = 0) with 

£(0,0) = -l+p|J_|-(l-p)J+ 

and ferrimagnetic minimum at n\ = p, 77.2 = or 771 = 
0, 772 = P with 



E(p,0) = - (1-p) 



P 



-p|J_|-(l-p)J+. 



Apparently, the ground state is the ferrimagnetic one 

(TOl = 1 — 2p, 7772 = 1 Or TOl = 1,7772 = 1 — 2p) if 

|J_|>l-p. 

The global (p, T) - phase diagrams presented in Fig. 
(JTJ for the whole exchange parameter space show that the 
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FIG. 1: The sections of the excange parameter space with different types of (p, T) phase diagrams shown in the inserts. Solid 
and dashed lines designate first- and second-order transitons correspondingly. The points' coordinates are: To — (0.658, 0.658), 
7\ - (0.619,0.619), T 2 - (0.418,0.690), T' 2 - (0.690,0.418), T 3 - (0.518,0.868), T3 - (0.868,0.518). 



FA phase can adjoin the paramagnetic one only at the 
point (p c , 1). Generally the transition into it is possible 
only from F and A phases and it can be either of second 
- or first - order type. 

The dotted curve in Fig. 0) divides the exchange 
parameter space in two regions with phase diagrams 
having four-phase point and only three-phase point 
at (p c , 1). This curve is defined by the equation 
3 tanh( J + ) tanh |J_ | = 1 and below it the stability re- 
gions of F and A phases always overlap resulting in the 
first-order transition between them. At J+ > 1 F-FA 
transition is always of the second-order type and the 
same is true for A-FA transition for |J_| > 1. When 
J + + |J_| < 1 the FA phase is always metastable. The 
regions 0-4 and l'-3' differ by the number of tricritcal 
points on the boundaries of FA phase at which the sec- 
ond - and first - order transitions meet. The lines be- 



tween these regions are given by cumbersome parametric 
equations, which we do not present here. 

Fig. I2| shows the typical temperature dependencies 
of the magnetic permeability, layers' magnetizations and 
heat capacity. They exhibit the characteristic mean-field 
anomalies at the second-order phase transitions. Note 
that magnetic permeability is not critical at the transi- 
tions into A phase and from F to FA phase, so it has a 
break in the first case and a jump in the second one. 
III. MAGNETIC PHASES IN EXTERNAL 
MAGNETIC FIELD AT T=0 

The study of phase states in magnetic field in the sys- 
tem considered is quite simple at T = 0. Here we can get 
analytical expressions for their stability regions. Indeed 
at T = we have from Eq. |J5J (putting again J = 1) 



F (m) = m 2 /2 — ( min 
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FIG. 2: Temperature dependencies of magnetic permeability, layers' magnetizations and heat capacity in the region 9 of Fig. 
©, (a)-p = 0.3, (b)-p = 0.7 



Here d (x) is Haeviside's step function. The form of 
the potential in Eq. JBJ implies that all its extrema are 
minima, so all solutions to eqations of state correspond to 
stable phases. Their stability regions and the equilibrium 
potential values are: 

1. F phases, m = ±1, I = 0, 

Hsign(m) > max(i7j, —1), 
H f = -1- J-, 

Fi = —1 — p J_ — (1 — p)J+ — 2Hsign(m). 



The stability regions of these phases and their magne- 
tizations m are shown in Fig. ijJJ for Hf > — 1 and in 
Fig. |0J for Hf < —1. Here the FA magnetization is 
shown for only one p value, while moving it inside the 
band bounded by dotted lines one can obtain its position 
for another p. Note that in Fig. |0J the cases with J± > 
are shown, which have the equilibrium F phase at all H, 
while A and FA phases appeare only as metastable ones. 

The first-order transitions can take place between co- 
existing phases at the fields in which their potentials be- 
come equal. Let us define the fields 



2. A phases, m = 0, I = ±1, 

\H\ < H a = 1 - J+, 

F 2 = -l+pJ_ + {l-p)J+ 



Hi =H a -1+p, 

H 2 =H f +p, 

H 3 = (l-p)H 1 +pH 2 . 



3. FA phases m = ±(1 — p), I = p, 

H sign{m) + \m\ > 0, 

H a < Hsign(m) + 2 |m| < Hf + 2, 

F 3 = -p 2 - (1 - pf + pj_ - (1 - p) J + - 2mH 



When Hi < < H 2 there are FA-FA transition at 
H = and FA-F transitions at H = ±H 2 while for < 
Hi < H 2 we have A-FA and FA-F transition at H = ±Hi 
and H = ±H 2 correspondingly. If H 3 < Hi just A-F 
transitions take place at ±H 3 . In other cases there are 
only F-F transition at H = 0, see Fig. (0J. The simple 
graphical procedure for obtaining these transition points 
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FIG. 3: Magnetic states in external field at T — for Hf > 
-1. (a) H a < 0, (b) < H a < 1, (c) 1 < H a . The dotted 
lines show the boundaries of the existence of FA states at 
varying p. Thin solid lines are used to determine the first- 
order transition points shown by the arrows. The tangents of 
the inclined lines are -1 and -1/2 except the line on Fig. 3(c) 
intersection of which with FA magnetization gives the A-F 
transition point. 



is shown in Fig. 

The graphs in Figs. ©, l@J give also a notion of the 
possible forms of hysteresis loops in slow varying peri- 
odic fields. In this case the transitions between magnetic 
states appear at the points where phases seize to exist 
since the period of field oscillations is usually much less 
than the decay time of metastable state. This implies 
also that hysteresis loops' forms could depend on the ini- 
tial state as well as the amplitude of the applied field. 
So for each graph in Figs.©, 0} we could have several 
loops. Yet the present static data are insufficient to pre- 
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FIG. 4: The same that in Fig. © for H f < -1. (a) H a < 0, 

-2 < H f < -1, (b) < H a < 1 

diet their exact forms in the cases when several stable 
states exist at a given field. Then one should consider 
kinetics as it is not generally evident to which state the 
system drops. For example, the exisence of small loops 
at low field amplitudes inside the large ferromagnetic one 
in Fig. Q can only be proved in the kinetic framework. 
Note also that the coercive field defining the boundary 
of ferromagnetic state is never greater than the one-layer 
coercive field being 1 in the present system of units. 

It is easy to envisage the evolution of the graphs in 
Figs. ©, (0J with the temperature rising. Thus the 
growth of permeability, see Fig. J2J , means that the lines 
depicting magnetizations become inclined while the field 
intervals where they exist become smaller. Finally the 
FA phases vanish at some temperature and above it only 
usual ferromagnetic and antiferromagnetic (double) loops 
would exist. 



IV. DISCUSSION 

The above theory can also be valid when instead of 
two monoatomic layers there are two slabs each hav- 
ing k such layers. Suppose the layers in the slabs 
are tightly bound by the nearest-neighbor ferromagnetic 
exchange (J') which is much stronger than intralayer (J) 
and random spacer (J) ones, J' ^ (kJ,J). Then all 
layers in the slabs would have parallel magnetizations at 
T <C J' so the present results hold after the substitution 
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J — > kJ . While such slabs are ubiquitous experimental 
objects, only the di ffere nt ones were usually studied in 
bilayers, see Refs. j ll'2y ;{0 llloOill l 1 .'ill 1 1| . Evidently this is 
due to expected triviality of results (just one magnetic 
state) for the identical slabs without randomness (yet 
see Ref. 0)- The present study show that random com- 
peting exchange can essentially complicate the variety of 
phase states in the simplest layered heterostructure. In 
case of different layers (slabs) the phase diagrams would 
have up to four ferrimagnetic phases as preliminary study 
shows. The variety of magnetic states in random multi- 
layer structures will be still more complex, and one may 
expect a fast phase alternation under small changes of 
external parameters. This could be seen in experiment 
as the region of spin-glass stated. 

Yet additional magnetic states in random heterostruc- 
tures could expand the possibilities of their technical ap- 
plications. Their study could also be fruitful for the gen- 
eral theory of disordered magnets. Thus the simplest 



system considered here allows one to conclude that the 
influence of competing random exchange does not come 
down to the changes of critical indexes. It can also cause 
the appearance of new magnetic phases and transitions 
between them under temperature as well as disorder vari- 
ations. It is natural to expect the existence of similar 
effects in the most studied objects with random compet- 
ing exchange such as spin-glasses with that distinction 
that they would have enormous multitude of phases and 
phase transitions. 
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